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We analyze renormalization of the K* band of n-doped epitaxial graphene on Ir(lll) induced by electron- 
phonon coupling. Our procedure of extracting the bare band relies on recursive self-consistent refining of the 
functional form of the bare band until the convergence. We demonstrate that the components of the self-energy, 
as well as the spectral intensity obtained from angle-resolved photoelectron spectroscopy (ARPES) show that 
the renormalization is due to the coupling to two distinct phonon excitations. From the velocity renormalization 
and an increase of the imaginary part of the self-energy we find the electron-phonon coupling constant to be 
~ 0.2, which is in fair agreement with a previous study of the same system, despite the notable difference in 
the width of spectroscopic curves. Our experimental results also suggest that potassium intercalated between 
graphene and Ir(l 11) does not introduce any additional increase of the quasiparticle scattering rate. 



I. INTRODUCTION 

Graphene is a fascinating two-dimensional material with 
linear electronic bands and linear dependence of the den- 
sity of states around the Fermi level. By field-effect or 
chemical doping (either n or p) Fermi level can be tuned 
to change graphene from a zero-gap semiconductor to a 
metal. " Intrinsic and doped graphene prove to be an inter- 
esting platform to study different aspects of many-body inter- 
actions in two-dimensional systems, either experimentally ~ 
or theoretically " ' . A particular accent is put on the electron- 
phonon coupling (EPC) in doped graphene. 

Adsorption and/or intercalation of alkali atoms in epitaxial 
graphenes can lead to a whole range of n dopings, including 
an extreme case where graphene' s van Hove point is brought 
to the Fermi level. Studies on KCg (Griineis et al. " ) and 
CaCg (Valla et al. ) report strong EPC anisotropy, with sig- 
nificantly smaller EPC constant A along the K-F, compared 
to the K-M direction. Despite the fact that these experi- 
ments were performed by intercalation of graphite, it had been 
demonstrated that the intercalation separates single layers of 
graphite in such a way that it shows all signatures of graphene. 
Some others have reported smaller value of A, " and no 
anisotropy in the EPC.'^^ Recent results also suggest that EPC 
progressively increases with doping. ^ 

Determination of the renormalization effects in «-doped n* 
band due to the EPC is complicated, in particular along the 
K-M direction, because of the non-linear behavior of the band 
close to the Fermi level. A simple procedure to establish the 
strength of electron-phonon coupling by a comparison of the 
quasiparticle velocity at the Fermi level and the velocity well 
beyond the phonon energy scale is shown to be compromised 
by a significant change of the quasiparticle velocity due to the 
nonlinearity of the n* band. Some other procedures rely on 
the bare-band dispersion.-" The choice of the bare-band may, 
however, influence an accurate determination of the magni- 
tude and shape of the self-energy. A self-consistent GW ap- 
proximation was applied in ab-initio density functional the- 
ory (DFT) calculation to model graphene's bare band that in- 
cludes all electron-electron correlations. In order to avoid 
any arbitrariness, determination of the bare band dispersion 
from the experimental data is desirable. Different approaches 



based on self-consistent procedures have been developed to 
determine the bare band, self-energy and ultimately, electron- 
phonon coupling strength. 

A self-consistent method has already been applied to 
the photoemission data of potassium doped graphene on 
Ir(lll). Several aspects of the low-energy quasiparticle dy- 
namics were addressed: renormalization of the K* band close 
to the Fermi level due to the coupling to phonons; phonon 
spectrum associated with the renormalization; the width of 
spectral lines in connection with the electron scattering rate. 
A model with a spectrum of five evenly spaced phonons par- 
ticipating in the coupling to graphene's n* band was used 
in the self-energy analysis. The EPC constant was found 
relatively small (0.28) and isotropic around K. Rather broad 
peaks were interpreted in terms of an increased electron scat- 
tering rate caused by the loss of translation symmetry induced 
by the incommensurability of the system graphene/Ir(lll). 
An ARPES study on another metallic system — graphene 
on a copper foil — questions all previous findings, propos- 
ing an order of magnitude smaller electron-phonon coupling 
strength. 

In this paper we present results of an ARPES study of 
highly ordered graphene on Ir(lll)-^^'-'* intercalated^ with 
potassium. We use maxima of momentum distribution curves 
(MDC) and energy distribution curves (EDC) to determine the 
exact dispersion of the n* band along the K-M direction. Peak 
positions and widths of MDCs are used in a self-consistent 
method to reconstruct the bare band Ef, and the coiTespond- 
ing 3xnL{E) and $HeE(£'). Both show, consistently with the 
spectral intensity A{E), that the renormalization is due to the 
coupling to two distinct phonon excitations. From the veloc- 
ity renormalization and an increase of 3mE(£') with energy 
electron-phonon coupling constant is determined. 



II. METHOD 

Photoemission spectra offer a wealth of information 
about many-body interactions in solids, in particular two- 
dimensional ones. The single-particle spectral function, that 
the intensities in the photoemission spectra are proportional 
to, is given by: 
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Here, Eh is the dispersion relation of the bare band, and E = 
IHcE + / 3mE its many -body correction, so-called self-energy. 
The latter is usually taken to depend only on energy, as its 
momentum dependence is considered weak. 

If a cut at a given energy E = E„, (momentum distribution 
curve) is made out of a two-dimensional map A(^,£'), nearly 
Lorentzian lineshape is obtained with a maximum at k„ such 
that 



E,„-Eb{k,„) - 5HeE(£',„) = 



(2) 



and a half-maximum at k,n — w^m and fc„, + WR,n (wim + w^,„ = 
2 w,„ is then full-width at half-maximum, FWHM), where 



3mT.{E,„) = Et,ik„,) -Eh{k,„ - wim) 



(3) 



as elaborated by Kordyuk et al.-''-'' Note that these relations 
do not rely on any specific dispersion of the bare band Ei,{k). 

The analysis of all the MDCs usually provides a set of value 
quadruplets (£,„, A:,,,, 2w,„, A„,) evenly spaced on £„,. If 
the lineshape of an MDC is not far from Lorentzian, one has 
WLm ~ WRm, and can set them both equal to w,„ — half-width 
at half-maximum, HWHM. The fact that each of the quadru- 
plets is supposed to satisfy equations (1), (2), (3) can help us 
determine the functional form of the bare-band Et,{k), and the 
two components of the self-energy 5HeE(£'), 3mE(£') in the 
range of momenta and energies covered by the spectrum. 

The procedure of extracting Ei,, 9^eE, 3mE from the experi- 
mental data assumes that (i) the components of the self-energy 
must, by causality, comply with the Kramers-Kronig rela- 
tion over the ±oo range of energies; (ii) the upper half of the 
spectral function not accessible to ARPES (E > fipermi = 0) 
is substituted by supposing particle-hole symmetry; (iii) the 
high-energy tails of the three functions are not affected by 
a finite energy window for which the ARPES spectrum is 
available.^ ''-^ 

If these criteria are met, the Kramers-Kronig transform of 
3mE to $HeE (and vice versa) 
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if considered as a convolution of two functions, is easily cal- 
culated going through the time domain by Fourier transforms 
(FFT, in the discrete case) 



[fHeE(£)]ppr = 



X pmE(£)]ppi. 



(5) 



FFT 



In some numerical packages this is already provided as a dis- 
crete Hilbert transform. 

The bare band and the components of the self-energy are 
usually reconstructed by assuming a polynomial (mostly lin- 
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Figure 1. (Color online) Self-consistent iterative algorithm that re- 
fines the bare band Ei, by enforcing the Kramers-Kronig relation be- 
tween 3vnL and SHeE. 



ear) form of the bare-band, from which DlzT.{E) is calculated 
by both eq. (2) and a Kramers-Kronig transform of 3mT.{E). 
Usually, not even the eq. (3) is used, but its expansion to the 
first order in w, 3m'E,{E,„) = h\'i,{Em) ■ w„,, where Vh(E) is the 
bare-band velocity. The two forms of 5HeE(£') are then com- 
pared, and the coefficients of the polynomial adjusted until the 
difference is minimal.-"'-''-^--** 

In this work we propose a simpler procedure, that avoids 
the need for a minimization of a functional, and achieves the 
self-consistency in just a few iterations. The procedure also 
appears to be more stable with respect to the noise present in 
the experimental data, and no denoising or smoothing by 
fitting to a function'^ is needed. 

The algorithm, shown in Fig. 1, starts by postulating a 
form, linear for example, of the bare band Et,{k). Forcing the 
function to pass through the experimentally determined point 
(kp, Ef) can help in faster convergence. Note that this con- 
straint is physically sound, as the renormalized and the bare- 
band are expected to intersect at the Fermi level. In the first 
step of the iterative procedure, 3mE is calculated by eq. (3). 
In the second step $HcE is obtained by a discrete Hilbert trans- 
form of 3mE. This is then by eq. (2) used to refine the bare 
band, first on a discrete set of points (step 3), then fitted to a 
polynomial or a function at wish (step 4). The iteration ends 
when there are no substantial changes in the functions (or pa- 
rameters) obtained. Further details on the procedure, the full 
code and exemplary input data are accessible as supplemental 
material in Ref /^J. 

Electron-phonon coupling strength can be extracted ei- 
ther from a steplike increase of the imaginary part of the self- 
energy, A3mE, as 



A = 



2 A3mE 
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or from the real part of the self energy, as 
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Although these two equations are strictly valid only at zero 
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temperature, they are especially good approximation in the 
case of graphene, where phonon excitations are high in en- 
ergy. 

Instead of calculating the derivative in eq. (7) directly, 
which is liable to errors due to the noise and scarcity of re- 
liable data in the very proximity of the Fermi level, one can 
make use of eq. (2), through A = —-^(Em^Ei,{k„,)) = 



A = ^-l (8) 

Here, v^, (vv) is the Fermi velocity of the bare (renormalized) 
band. 



III. EXPERIMENTAL 

The experiments were performed at the APE beamline at 
ELETTRA. Iridium single crystal of 99.99% purity and sur- 
face orientation better than 0.1° was used. The substrate was 
cleaned by several cycles of sputtering with 1.5 keV Ar"*" ions 
at room temperature followed by annealing at 1600 K. Clean- 
ness and quality of Ir(lll) were checked by ARPES (exis- 
tence and quality of iridium surface states) and low energy 
electron diffraction (LEED). The base pressure was 510"^ Pa. 

During the growth of graphene the surface was exposed to 2 
Langmuir (L, 1L=1.3-10^ Pa-s) of ethene at 300 K and subse- 
quently heated to 1470 K. In order to ensure perfectly oriented 
graphene and the coverage of the whole surface, the proce- 
dure was repeated up to 12 times, finishing with simultaneous 
exposure to ethene and heating. ' Potassium was added at 
room temperature by evaporation from a getter source, until a 
clear 2x2 structure emerged in LEED and no additional dop- 
ing of graphene bands could be achieved. 

ARPES spectra have been collected by a Scienta SES 2002 
analyzer. For this experiment we used photon energy of 
40.5 eV, in s and p polarization of the field. Typical spot size 
on the sample was 50x100 jim. The energy resolution of this 
setup was about 12 meV, and the angular resolution 0.1°. The 
azimuths were checked according to the orientation of spots 
in LEED and the emergence of the Dirac cone and replica 
bands in ARPES for the plain graphene on Ir(l 1 1). During 
the measurements the temperature of the sample was held at 
80 K. 



IV. RESULTS 

Figure 2 shows photoemission spectrum of potassium- 
intercalated graphene on Ir(l 1 1) around the K point along the 
F-K-M direction. The spectrum shows a discernible asymme- 
try of the intensity, such that the part of the spectrum of the n* 
band along K-M (k > 1.7 A"^) exhibits much stronger inten- 
sity than the part dispersing along K-F {k < 1.7 A"'). Notice 
that iridium surface state at the Fermi level, also present in 
the spectrum for graphene/Ir(l 1 1), is not quenched by the 
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Figure 2. (Color online) ARPES spectrum (excitation energy 
40.5 eV, s polarization, scan direction F-K-M) showing the intensity 
of photoemission from the Dirac cone in graphene/Ir(lll) saturated 
with potassium. 



intercalation of potassium. Potassium doping, however, has 
rather strong effect on the K bands of graphene — shifts the 
Dirac point to higher binding energies and renormalizes the 
band dispersion just below the Fermi level. 

The position of the Dirac point, being defined as a single 
point in momentum space where the n bands cross each other, 
is not straightforward to determine. As can be seen from 
Fig. 2 there is no such well-defined crossing point here. This 
is presumably due to a band gap, as a study by Varykhalov 
et al. has shown that some metallic dopands do induce a 
band gap at the Dirac point. We estimate the width of the band 
gap to be 0.3 e V with the position of the Dirac point at 1 .35 e V 
below the Fermi level, which is, within a tenth of an eV, equal 
to the values previously obtained for graphene/Ir(l 1 1) and 
some other systems, such as graphene/SiC and graphite . 

The shift of the Dirac point to higher binding energies 
causes an increase of the Fermi surface, which is for higher 
doping levels characterized by trigonal warping; i. e. an effect 
when a transformation of the constant energy maps from cir- 
cular to trigonal shape takes place. ■ The trigonal warping 
is clearly associated with different values of kp for the n* band 
dispersing from K to M or F. We have determined kf along 
the directions K-F and K-M to be 0.18 A"^ and 0.27 A"' re- 
spectively, and used these values to estimate the doping level 
of about 0.05 electrons per unit cell. 

A prominent feature of the spectrum in Fig. 2 is the renor- 
malization of the K-M branch of the n* band just below the 
Fermi level. No such renormalization is as obvious for the 
branch K-F. This kink-like change of the band dispersion is 
detailed in Fig. 3. 

As we have pointed out, two parts of the spectrum in Fig. 
2 show a considerable difference of spectral intensity. The 
change of the light polarization can additionally alter this in- 
tensity ratio. We used p polarization of the incident light in 
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Figure 3. (Color online) ARPES spectrum of potassium intercalated 
graphene/Ir(lll) (excitation energy 40.5 eY, p polarization, scan di- 
rection F-K-M) along with a few selected MDC and EDC cuts at 
given energies and momenta. 



order to extinguish spectral features along K-F entirely, leav- 
ing only the part of the n* band along the K-M direction visi- 
ble. The origin of the diminishing spectral intensity has been 
explained by Gierz et al. The change of the light polariza- 
tion strongly enhances signal to noise ratio which shows to 
be essential in detailed spectrum analysis, as will be demon- 
strated in the following. Simple inspection of Fig. 3 places 
the phonon induced dispersion kink at around 200 meV be- 
low the Fermi level. Notice how the dispersion kink is now 
even more pronounced, accompanied by a strong drop of the 
spectral intensity at the kink. In the following we present a de- 
tailed analysis of the spectrum along the K-M direction, and 
only state the results for the K-F direction for which the anal- 
ysis is presented in the supplemental material, Ref. 29. 



V. ANALYSIS AND DISCUSSION 

A set of MDC and EDC cuts in Fig. 3 illustrates the kind of 
analysis made for each and every slice of the ARPES spectrum 
shown. MDCs are generally characterized by a continuous 
change of the peak position without any substantial change of 
the lineshape, apart for the width of the peaks. However, in 
the narrow region around the kink, the EDCs clearly exhibit 
double peak structure, that can be fitted with two Lorentzians. 
The energy splitting between the peaks is found to be around 
60 meV. Farther away, +70 meV from the kink, EDCs can 
be described by a single Lorentzian. A small peak visible 
in EDCs just below the Fermi level is associated with the S i 
surface state of Ir(ll 1). 
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Figure 4. (Color online) Peak positions extracted from MDC cuts 
(circles) and EDC cuts (squares) of the spectrum shown in Fig. 
3. Thick line shows the bare-band obtained from the self- 

consistent procedure. 



Note that the full width at half maximum of the MDC at an 
energy just below the Fermi level (Fig. 2) is 0.022 A"', which 
is even slightly smaller than the width measured for bare 
graphene on fr(lll)"' and comparable to non-intercalated 
graphene on SiC^. The measured width of 0.022 A"' is 
substantially smaller than the previously reported value of 
0.095 A"' for the same system. This questions the conclu- 
sion by Bianchi et al. that the doping of graphene/Ir( 111 ) by 
potassium increases the electron scattering rate due to the loss 
of translation symmetry induced by the incommensurability 
of graphene and Ir(lll). Our data imply that the intercala- 
tion of potassium, observed as a disappearance of graphene's 
moire superstructure, does not increase the electron scattering 
rate compared to bare graphene on Ir(l 1 1). Therefore, we can 
conclude that potassium intercalated into graphene/Ir(l 1 1) 
system does not act as an additional scattering center. As we 
shall demonstrate later, the measured value of the MDC width 
translates into rather big quasiparticle scattering time. The 
MDC width increase (to 0.042 A"') measured below the dis- 
persion kink (Fig. 2) fits into the picture of phonon-induced 
renormalization of the n* band. 

Figure 4 summarizes the peak positions obtained from 
MDC and EDC sets of spectra from Fig. 3. Open circles 
represent peak positions obtained from the analysis of MDCs, 
while squares are obtained from EDCs. The fact that the two 
dispersions do not coincide, especially in the region above the 
kink is easy to understand when one takes into account that 
the spectral intensity of the band sharply changes, and that 
what is a local maximum k,„ in a horizontal cut along E = Em 
does not have to be a local maximum in a vertical cut along 
k = k,„- The data below and above the kink are fitted with two 
linear functions with noticeably different slopes. The data far 
away from the Fermi level can be fitted with a linear curve 
and the corTesponding band velocity equals to 0.7 • 10^ m/s. 
This represents a reduction of about 30% compared to the n 
band of nearly neutral graphene on Ir(lll), and is consistent 
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with the predicted reduction of the Fermi velocity, as calcu- 
lated by Park et al. The set of MDC data between the Fermi 
level and the kink is also fitted with a linear function defined 
by the band velocity (which is in this case the Fermi velocity 
as well) Vr{Ep) — 0.46 • 10^ m/s. Note that this value is only 
about half the Fermi velocity of the n band in bare graphene 
on Ir(lll) and the n* band on SiC^. This renormalization 
of the velocity can lead to an overestimate of the coupling 
strength to phonons for the simple reason that the change of 
the electron velocity near Ef is largely due to the curvature of 
the n* band induced by other interactions and only to small 
extent by electron-phonon coupling. 

In order to take into account the nonlinearity of the n* 
band, the bare-band function Ei, (k) was self-consistently re- 
constructed from the MDC data, by the recursive procedure 
described above. The n* band itself has been measured 
far below the phonon-induced kink, where the contribution 
to IHeE can be neglected and the bare-band approaches the 
measured. This gives confidence to the functions obtained, 
as the problem of tails is greatly avoided. The corrections 
of Ei„{k) leading to Ei,{k) being mostly due to low -energy 
phonon excitations, the function E[,{k) can be considered as 
one that takes into account all electron interactions except 
those with phonons. 

The derivative of the bare function at the Fermi level gives 
the Fermi velocity for the bare band Vh{Ef) = 0.56 • 10^ m/s. 
The obtained values of band velocities, at the Fermi level and 
far away from the kink, are in very good agreement with the 
values of the band velocities for highly doped graphene, re- 
ported by Siegel et al. 

To summarize, it is clear that a large portion of the velocity 
reduction as we approach the Fermi level comes from the non- 
linear nature of the bare n* band. The modification of the 
Fermi velocity due to the coupling to phonons is accordingly 
rather small, around 0. 1-IO'' m/s, which is only a 10% change. 

The ratio between the Fermi velocity of the bare band 
and the electron-phonon-induced renormalized velocity at the 
same energy gives, by eq. (8), the value of the electron- 
phonon coupling strength. Accordingly, the value of X along 
the K-M direction equals to 0.22. The same analysis on the 
K-F direction gives A = 0.19. ' Previously, Bianchi et al. 
obtained A =0.28 for both directions. 

We use MDCs extracted from the spectrum shown in Fig. 
3 to plot (a) imaginary part of the self-energy, 3mE(£') (b) 
real part of the self-energy, D\eT.{E) and (c) spectral inten- 
sity of the MDC peaks, A{E). Figure 5 shows the results. 
The functions plotted in black were calculated from formulae 
(l)-(3) using the bare-band from our self-consistent iterative 
procedure. The functions shown in red illustrate the degree 
of self-consistency achieved; the one for ^ttiEkk in Fig. 5(a) 
was obtained by Kramers-Kronig transformation of the data 
for *KeE from Fig. 5(b), and equally, the one for fHeEKK from 
the data for 3mE. The overall agreement is quite good. We 
find the data for FWHM and consequently for 3mE more re- 
liable, the noise being pretty low and the relative error being 
smaller in determining the width of the MDCs than the peak 
positions, from which $HeE is calculated. The effects of ex- 
perimental smearing also have a greater impact on the fine 



details of MDC peak positions, adding only a constant back- 
ground to FWHMs and 3mE. We suspect that the low-energy 
shoulder in 5HeE (and the associated step in 3mLKK) are for 
these reasons less pronounced than the corresponding ones in 
D^cLkk (and JmE). The difference in experimental and calcu- 
lated spectral intensities in Fig. 5(c) at the high-energy side 
probably comes from the fact that the calculated renormalized 
band Ei,{k,„) +£HeE(£'„,) does not come close enough to the 
measured band E,„{k„,), as the increase of the spectral func- 
tion would come from one approaching the other (see Eq. 1). 
We suspect this could be a still-observable consequence of the 
missing high-energy tails. 

Generally, a step-like increase of 3mE(£') and a maximum 
of lKeE(£') at the same energy is associated with electron cou- 
pling to some excitation, in this case phonons. As pointed out 
by McChesney et al. the coupling to a phonon should have 
a profound impact on the energy dependence of the spectral 
intensity A (£■) as well. 

Instead of one, all displayed spectral parameters show con- 
sistently two features that can be accordingly associated with 
the coupling to two phonons, one at 170 meV and the second 
one at 75 meV. Arrows in Figs. 5 (a)-(c) indicate the energies 
of these two phonons. 

The 3mE(£') (Fig. 5a) shows a distinct step-like increase 
between 0.05 and 0.1 eV below the Fermi level with the initial 
value of 39 me V which reaches a local maximum of 52 me V at 
0. 1 eV. This increase of 3vnL{E) is associated with the phonon 
of energy 75 meV. Further increase of 3mE(£') up to 85 meV 
induced by the coupling to a phonon of 170 meV takes place 
between 0.13 and 0.21 eV. Below that, '3xnL{E) stays constant 
up to 0.4 eV, but then starts to increase with energy, due to 
other contributions to the lifetime (electron-electron, electron- 
plasmon). 

Electron-phonon coupling constant can be extracted from 
the change of 3xnL{E) by eq. (6). ■ The total increase of 
3mE(£') equal to 46 meV, translates to A = 0.17 ±0.01 along 
the K-M direction. A similar analysis of 3xnL{E) along the 
K-F direction gives A = 0. 18 ± 0.02. 

5HeE(£'), consistently with 3mE(£') andA(£'), shows apeak 
at 170 meV and a distinct shoulder at a lower binding -energy 
side, which clearly implies an existence of two excitations that 
contribute to the renormalization of the bare band dispersion. 

Using the imaginary part of the self-energy in the proxim- 
ity of the Fermi level, we determine the photohole scattering 
time, T = 2 ymL ' bigger than 10 fs. This is the same 

value as obtained for high quality graphene on SiC. Note 
that this value is not to be compared to the scattering time 
in mobility measurements (^ 350 fs), as those include some 
more scattering mechanisms. 

As pointed out, the energy dependence of the spectral inten- 
sity A{E) also indicates the existence of two distinct phonon 
excitations that couple to electrons. McChesney et al. have 
demonstrated a high sensitivity of A{E) to many-body inter- 
actions, drawing attention to the fact A{E) can access even 
a faint contributions to the self-energy with a sensitivity even 
better than 3mE(£') can provide. In agreement with this model 
we observe two shoulders in A{E) where an onset of each 
shoulder corresponds to a phonon excitation (see Fig. 2 in 
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Figure 5. (Color online) A set of self-consistent functions obtained from MDC cuts of the spectrum shown in Fig. 3 and the bare-band shown 
in Fig. 4. (a) black line: DimZ calculated by eq. (3); red line: 3ml,Kx as a Kramers-Kronig transform of 9^cS; (b) black line: IHeZ calculated 
by eq. (2); red line: ^e'LfcK as a Kramers-Kronig transform of JmZ. (c) spectral intensity A{E) from the experiment (blue line) and calculated 
from the self-consistent data by eq. (1) (black line). Arrows indicate the energies of phonons that induce band renormalization. 



Ref. 28). Summarizing the features that have been observed 
in three parameters: 3m.'L{E), 9\eL{E) and A(£'), we can con- 
clude, with a high degree of consistency, that the phonon ener- 
gies that induce the renormalization of the dispersion of the n* 
band of doped graphene along the K-M direction have ener- 
gies equal to 170 meV and 75 meV. These two phonons can be 
associated with an optical phonon (transverse or longitudinal) 
around K and an acoustic phonon, respectively. According 
to Gonzalez and Perfetto , in non-interacting graphene these 
phonons should be in-plane oscillations as the symmetry does 
not allow coupling of the out of plane phonons to electrons in 
graphene. However the presence of the substrate might break 
the 2D symmetry of graphene and hence allow the coupling 
of both in-plane and out-of-plane oscillations to electrons in 
graphene. 

Some theoretical calculations also support the notion 
of a two-phonon-spectrum that induces renormalization of 
graphene bands close to the Fermi level. According to Calan- 
dra and Mauri , the 7t* band of doped graphene should show 
two kinks in the dispersion (accordingly accompanied by two 
steps in the MDC linewidth), one at 195 meV being attributed 
to Ai mode and another at 160 meV which corresponds to 
a twofold degenerate E2g mode. Previous measurements of 
Bianchi et al. with graphene on Ir(lll) showed signatures 
in the band dispersion and widths of spectral curves of appar- 
ently only one phonon. However, their self-consistent anal- 
ysis showed that the experimental data can be modelled by 
contributions of five Einstein oscillators with energies that are 
evenly distributed over the range from 21 meV to 190 meV. 

Interestingly, the pattern of coupling similar to the one we 
observed in graphene/Ir( 11 1) is also reported for CaCg where 
the contributions to the self-energy come from two phonons, 
one at 75 meV and the other at around 160 meV. This agree- 
ment supports the notion that intercalated atoms (K or Ca) do 
not participate in the coupling of the n* band to phonon modes 
apart from the doping effect. 

The energy splitting between the EDC peaks has already 
been observed in materials with strongly coupled electrons 



and phonons. ' Although the coupling in graphene is not 
as large, its two-dimensionality is, according to Badalyan and 
Peeters , a possible reason for an enhancement of the ef- 
fective coupling in the vicinity of the phonon energy. An 
electron-phonon complex quasiparticle forms, giving rise to 
a strong modification of the Dirac spectrum, seen as the 
branches of the n* band in the neighborhood of the phonon 
emission threshold. 

The value of the electron-phonon coupling constant ob- 
tained in this work (^ 0.2) is somewhat smaller than the one 
obtained by Bianchi et al. for the same system (^ 0.3). 
Nevertheless, the value derived from our spectra does not sup- 
port the notion of Siegel et al. that the coupling in n-doped 
graphene on a metallic surface should be still a few times 
smaller (^ 0.04). Given sharp enough spectra, low in back- 
ground intensity and noise, a pronounced kink in the n* band 
dispersion, and especially a steplike increase of the spectral 
width at the phonon energy, must alone be a clear sign of a 
non-negligible strength of interaction, whatever the choice for 
the bare band is. 



VI. CONCLUSIONS 

In conclusion, we have analyzed various parameters from 
high resolution ARPES measurements — peak positions and 
widths of MDCs, imaginary and real part of the self-energy 
£(£■), spectral intensity A{E) — all consistent with the no- 
tion that two phonons contribute to the renormalization of the 
band dispersion in graphene on Ir(ll 1): a high energy phonon 
at 170 meV and a low energy phonon at 75 meV. The cou- 
pling constant associated with these two phonons is around 
0.2, which is similar to a previous finding for the same sys- 
tem. We have also found that the intercalation of potassium 
up to saturation does not increase the scattering rate of the 
photohole in the n* band. Due to the perfectness of the KJ- 
graphene/Ir( 111 ) structure the quasiparticle (photohole) scat- 
tering time can be posed exceptionally high, above 10 fs. 
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SUPPLEMENTAL MATERIAL 
Formulae 

Under the assumption that the self-energy does not depend on momentum, equations (2) and (3) from 
our article are exact whatever the bare -band function Eh{k) is. However, in order to use the linearized 
version of eq. (3) — 3mL{E^) = hvb{E^) ■ w„„ where vt{E^) = j^^Ehik) \k=k„,{E„,) — the measured band 
should be sufficiently sharp (wm sufficiently small), and the bare band piecewise linear (giving a symmetric 
peak, wiin = WRm = Wm)- Even in our case of very narrow MDC peaks (0.010 A'^<w^<0.022 A"^, or 
0.007 A"^<Wot<0.018 A"^ when the experimental resolution is deconvolved) this is only partially fulfilled. 
The figure below illustrates the relations for JmE, 9^eE and E^, in the case of real experimental data on 
E,n{km) and FWHM (2w„,). Also shown is an MDC curve for = -0.13 eV with FWHM of 0.026 A'^. 

In the following, we shall not use the linearized version of eq. (3), in spite of rather narrow peaks, and 
only make use of wtm = WRm = Wm, as our Lorentzians are quite symmetric. 




(k-kp)lA 



10 



Method 



Instead of finding Ei,{k) by cleverly trying 
the parameters that define the function in a 
method commonly used before (Method I on 
the right), using Levenberg-Marquardt algo- 
rithm, for example, the idea is to employ 
a loop algorithm that happens to be self- 
correcting and relatively easy to implement 
(Method II on the right). 
Here, D^eE(£'„,) obtained from the Kramers- 
Kronig transform of IJmE is used to calcu- 
late Eh from the measured renormalized dis- 
persion Em{k,n) on a discrete set of points. 
Then, by least-squares fitting, a functional 
form Ei,{k) is acquired, and used to calculate 
3mL in the next iteration. 



METHOD I 



r 



ansatz - 



9^eS(£;) =E- Eb{h 



3mi:{E) = hvb{E)w 

jKramers-Kronig 

minimize the functional on | 

<KeE(£;) - ^zY.kk{E) overall E's 
by (cleverly) trying {man y) different Eb{k) 



METHOD II 



ansatz Eh{k) 



: through [kF.Ep) 



3mi:{Em) = Eb{k 

I Kramers-Kronig 



"1 



Eb{km) = Em{km) - ^^c'Sk K (Emikm)) 
I kj^, — >■ k polynomial fit, through {kF^Ep) 

Eb{k) 
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Convergence 

The figure below shows the convergence of the bare-band function. The iteration starts by postulating a 
linear band (io). Already the second iteration (i2) gives the band indistinguishable from those obtained in 
the following iterations. Also shown is how the coefficients of a third-order polynomial for Ei, 

Ebik) = + cik + C2k^ + C3k^ 

converge oscillatory through the course of iterations. 

Notice that the renormalized and the bare-band merge quite close to the kink, giving confidence that the 
influence of the missing tail is negligible. 
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Code 

Two procedures to be used in Wavemetrics IgorPro data analysis software are presented. The first pro- 
cedure, scKK^setup, is used for initialization, with properly scaled one-dimensional waves containing data 
on the positions (km) and widths (fwhm) of MDC Lorentzian curves; examples are given in the next section. 
The second procedure, scKK^ite ration, is called for several times, until the convergence in the output waves 
(scEb, scEb_f it, scims, scReS. . . ) is achieved. 

We find it convenient to put (kf,EF)=(0,0), as that point can be read-off the spectrum with high enough 
accuracy, and make the bare-band function pass through it. 

Instead of a polynomial, any fitting function (e.g. a tight-binding formula) should be easy to implement. 
As we are usually not interested in Efy, but in the results obtained for 3xxiL and 5^eE, polynomial fitting is 
sufficient. 



// 

// Author: Ivo Pletikosic, ivo . pletikosicOif s . hr, November 2011 
// 

// sc = self -consistent; KK=Kramers-Kronig 



// first, create waves needed for the iteration procedure 

// wave kni(Em) contains the experimental dispersion, i.e. maxima of the spectral function 
// wave fwhm(Em) contains widths of Lorentzian peaks at a given energy 
// X axes should be properly scaled! 
// energy is measured from the Fermi level, (E-EF) 
Function scKK_setup( km, fwhm) 
Wave km, fwhm 

// (km-kF) as a function of Em 
Duplicate/0 km,scKKm 

// putting the origin of the k axis to kP (a point where Em=0) 
Variable kF=km(0); scKKm=km-kF 
Duplicate/0 fwhm.hwhm; hwhm*=0.5 

// sample ReS at the same energies as KKm 
Duplicate/0 scKKm,scRe5 

// sample ImS at the same energies as hwhm 
Duplicate/0 hwhm.scImS 

// discrete bare band Eb; Em(km); and an auxiliary wave 
Make/0/N=(numpnts{scKKm) ) scEb,scEm,Em_scaling 

// scaling from the first/last k-point of the experimental dispersion 
SetScale/I x,scKKm[0] ,scKKm[numpnts{scKKni) ] , "A\S-1\H" , scEb, scEm 
SetScale y,0,0, "eV" , scEb, scEm 

// this will be used in the interpolation to get Em{km) from km(Em) 
SetScale/I x,leftx(scKKm) ,leftx(scKKm)+deltax(scKKm)* (numpnts(scKKm) -1) , Em_scaling 
Em_scaling=x 

// Em as a function of (km-kF) 
scEm=interp(x, scKKm,Em_scaling) 

// zeroth approximation to Eb: 

// a line through (kF,EF)=(0,0) and the lowest point of the experim. dispersion 
Make/0 cw={0,leftx(scKKm}/scKKm[0]} 
scEb=poly(cw,x) 

// this will become symmetrized and resampled ImS 
Make/0/N=(4+numpnts{scImS) ) symmlmS 
SetScale/I x,leftx(scImS) , -leftx(scImS) , "eV" , symmlmS 

// this will become symmetrized and resampled ReS 
Make/0/N=(4*numpnts(scReS) ) symmReS 
SetScale/I x,leftx(scReS) , -leftx(scReS) , "eV",symmReS 
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// iteration function 

// call it several times until there's no substantial change in the output waves 
// polyN is the degree of the band fitting polynomial [2, 3, 7] 

Function scKK_iteration(polyN) 
Variable polyN 

Wave scKKm, scEb, scReS, scImS 

Wave Em, hwhm, cw, symmReS, symmlmS // some auxiliary waves 

PauseUpdate; 

// start from ImS 
scIm5=poly(cw,scKKm(x) ) - poly(cw, scKKm(x) -hwhm(x) ) 
// or 

//scImS=poly(cw, scKKm(x)+hwhm(x) ) - poly ( cw, scKKm(x) ) 
// or 

//scltn5=0.5*(poly{ cw, scKKm{x)+hwhm(x) ) - poly(cw,scKKm(x) -hwhm(x) ) ) 



// K-K transform ImS to HReS 



// symmetrize 
symmImS=(x<=Q) ? scImS(x) : scImS(-x) 

// ReS from the K-K transform of ImS 
HilbertTransform/DEST=scHReS symmlmS 
scHReS*=-l 

SetScale/P x, leftx( symmlmS) ,deltax( symmlmS) , "eV" ,scHReS 



// ansatz for Eb in the next iteration 



// bare band as a function of km (therefore, of all the discrete k's present in the spectrum) 
scEb=Em(x)-scHReS{Em{x)) 

// fit Eb to a polynomial {could be any other function), as we'll need some extrapolation 



// initial values of the fitting polynomial coefficients: all 
Make/0/N=(polyN+l) cw=0 

// keep c0=0, i.e. make the polynomial zero at (kF,EF)-(0,0) 
String fmtStr; sprintf fmtStr, "l%0*.Qf " polyN, 0; 

// functional (i.e. smooth) form of the bare band 
Duplicate/0 scEb, scEb_fit 

// the output of this cmd will be used for the next iteration 
CurveFit/Q/H=fmtstr/NTHR=l poly (polyN+1), kwCWave=cw, scEb /D=scEb_fit 

// print fitting coefficients 
print cw 



// and just to check things out 

// ReS from the very definition 

// poly(...) represents a smooth function for Eb from the previous iteration 
scReS-x - poly(cw,scKKm(x) ) 

// antisymmet rize 
symmReS={x<=0) ? scReS(x) : (-scReS(-x)) 

// ImS from the K-K transform of ReS 
HilbertTransform/DEST=scHImS symmReS 
SetScale/P x, leftx( symmReS) ,deltax( symmReS) , "eV" ,scHImS 

// offset, as HilbertTransform seems to make mean=0 
scHImS+=mean(scImS) 
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Data 



Examples of input waves for the self-consistent iterative procedure, as extracted from MDCs of an 
ARPES spectrum analyzed in our paper. 



IgorPro file km . itx 

(here downsampled to A = 0.012 eV from the original A = 0.005 eV) 



IGOR 




















WAVES 


km 


















BEGIN 




















I 


797 


I 


799 


1 


801 


1 


804 


1 


806 


I 


809 


I 


812 


1 


814 


1 


817 


1 


819 


I 


822 


I 


824 


1 


827 


1 


83S 


1 


834 


I 


836 


I 


838 


1 


841 


1 


844 


1 


846 


I 


849 


I 


852 


1 


854 


1 


857 


1 


859 


1 


862 


I 


864 


1 


867 


1 


87S 


1 


872 


1 


875 


I 


878 


1 


88S 


1 


883 


1 


886 


1 


889 


I 


891 


1 


893 


1 


895 


1 


898 


1 


900 


I 


903 


1 


905 


1 


907 


1 


910 


1 


913 


I 


918 


1 


922 


1 


927 


1 


932 


1 


936 


I 


94S 


1 


944 


1 


947 


1 


951 


1 


955 


1 


959 


1 


963 


1 


968 


1 


971 


I 


975 



















END 

X SetScale/P x 
X SetScale y S, 



-0.7,0.012, "( 
S,"l/A", km 




IgorPro file f whm .itx 

(here downsampled to A = 0.012 eV from the original A = 0.005 eV) 



IGOR 




















WAVES 


fwhm 


















BEGIN 




















S 


B464 





0458 


B 


B453 





0434 





B430 


G 


B445 





0435 


B 


B442 





0423 





B414 


B 


B420 





0407 


B 


B418 





0419 





B415 


B 


B405 





0414 


B 


B420 





0416 





B4B8 


B 


B398 





0390 


B 


B4B0 





0405 





B4B8 


B 


B397 





0398 


B 


B399 





0403 





B4B5 


B 


B408 





0408 


B 


B413 





0406 





B399 


B 


B419 





0421 


B 


B411 





0415 





B418 


B 


B433 





0436 


B 


B428 





0382 





B361 


B 


B315 





0290 


B 


B282 





0273 





B279 


B 


B278 





0262 





0243 





0233 





B230 


B 


B228 





0224 





0245 





0247 





B246 


B 


B222 


















END 




















X SetScale/P 
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A few questions to answer 

• Could the K-F direction be treated with the same rigor as the K-M direction ? 

In the experimental geometry we had, the K-F direction (fc||<1.7 A"^) is not visible in high-intensity spectra 
recorded with the p polarization of light. This is illustrated in the following two raw images taken with 
the same excitation energy (40.5 eV) and in the same energy and emission-angle window. One can easily 
notice a high noise level in the s polarization spectrum where a part of the K-F direction is seen. The 
photoemission intensity is about 30 times lower in the s polarization spectrum. 




kn/A"'' k||/A" 
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• So, what are the results for the K-F direction? 

We shall analyze the spectrum recorded in the s polarization of the excitation light, shown below. The 
whole K-F branch of the doped Tt* band is visible, continuing with a part of the K-M branch. A surface 
state of Ir(l 11), touching the renormalized part of the n* band is also visible. It influences the analysis of 
the MDC peak positions and widths in a small range of energies, between -25 meV and -70 meV. These 
data will have to be interpolated. 




We now show the data on FWHM and the position of MDC peaks extracted from the above spectrum, 
along the K-F branch of the band. The widths are comparable to those on the K-M branch, but the noise 
in the data is obviously larger: 




Applying our self-consistent recursive procedure to the data from the K-F branch, we obtain the bare 
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band (shown above, along with the dispersion of the renormalized band), and the parts of the self energy, 
E = 9ieE + /amE: 
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• Can you simulate a spectral function with the same lambda (X = 0.2j in the K-M and K-T direction. 
It would be interesting to see if the different dispersion in the two directions is responsible for the kink 
being more obvious in the KM direction. 

We simulate this by first recognizing that the same e-ph coupling parameter A = 0.2 for the phonon of 
energy (Oph = 170 meV means, by equation (6) of our article, the same step in the imaginary part of the self 
energy. This would be AUmE = 55meV for both branches. The real part of the self energy is obtained by 
the Kramers-Kronig transformation of the imaginary part, assuming particle-hole symmetry: 
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We then take two linear bare bands, one mimicking the dispersion of the K-F branch (band velocity 
1 • 10^ m/s), the other the dispersion of the K-M branch (band velocity 0.6 ■ 10^ m/s), and simulate the 
spectral function by equation (1) from our article: 
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simulation 
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The renormalized bands can be either extracted from the spectral function following the dispersion of 
MDC peaks, or calculated by equation (2). This is done self-consistently, starting with E^^ = Ei,{k) + 
y{zL{Eh{k)) and iterating over E]:^^ = Et,{k) + 9^eE(£'^). The result can be immediately compared to the 
experimental data. Even the simple linear band closely follows the measured band in a wide range below 
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the kink, as well as its renormalized part follows the band above the kink. And indeed, the kink is more 
obvious for the slowly dispersing band, i.e. the band simulating the dispersion along the K-M direction. 
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